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Abstract 

A theory of correlations between N photons of given frequencies and detected at given time delays 
is presented. These correlation functions are usually too cumbersome to be computed explicitly. 
We show that they are obtained exactly through intensity correlations between two-level sensors 
in the limit of their vanishing coupling to the system. This allows the computation of correlation 
functions hitherto unreachable. The uncertainties in time and frequency of the detection, which 
,__! are necessary variables to describe the system, are intrinsic to the theory. We illustrate the power 

o 

\Q of our formalism with the example of the Jaynes-Cummings model, by showing how higher order 

photon correlations can bring new insights into the dynamics of open quantum systems. 
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Photons emerged as a theoretical concept to explain fundamental properties of the elec- 
tromagnetic field, such as the relationship between the energy of light and its frequency, 
thermal equilibrium of light and matter or the photo-electric effect. With the advances 
in the generation, emission, transmission and detection of photons, quantum systems are 
increasingly addressed at the single photon level and there is a pressing need for general- 
izations as well as refinements of the theory of photo-detection [1]. For instance, photon 
correlations combining both their frequency and time information are now routinely mea- 
sured in the laboratory. These experiments have proven extremely powerful in characterising 
quantum systems such as a resonantly driven emitter [2-4], the strong coupling of light and 
matter [5-7], to perform quantum state tomography [8], to monitor heralded single photon 
sources [9] or to access spectral diffusion of single emitters [10]. 

At this level of fine control of the attributes of the quantum particles, one needs a the- 
oretical description significantly more involved than general mathematical statements, such 
as the Wiener-Khinchin theorem which assumes abstract and unphysical properties of the 
light field. Eberly and Wodkiewicz, for instance, have shown how the physics of the de- 
tector needs to be included if a more realistic description of the light field is required [11]. 
In general, the more detailed is the characterization of a quantum system, the more neces- 
sary it becomes to describe its measurement. A bridge between the quantum system and the 
observer can be made with the so-called input-output formalism: the photons inside the sys- 
tem, say with operator a (we consider a single mode for simplicity), are weakly coupled to an 
outside continuum of modes, with operators A u (corresponding to their frequency u). In the 
Heisenberg picture, the output field allows to compute the time-dependent power spectrum 
of emission as the density of output photons with frequency oj\ at time 7\, i.e., S^(u>i, Ti) = 
(A' (Ti)A ul (Ti)) . This quantity is physical only if the uncertainties of detection in both 
time and frequency are jointly taken into account [11]. Mathematically, this amounts 
to adding two exponential decays in the Fourier transform of the time-autocorrelation 
4'WTi) = ^//^ o «e-^( Tl - t i)e-^( Tl - i 4)e^(*4-ti)( a t(t' i ) a (^)) where I\ is inter- 
preted as the linewidth of the detector. This so-called physical spectrum reduces to the 
Wiener-Khinchin theorem in the steady state and in the limit r x — > 0. 

Extending this result for the detection of two photons was initially motivated by the 
Aspect et al. experiment [2] of resonance fluorescence in the Mollow triplet regime [12], 
where the peaks of the triplet were found to exhibit strong intensity correlations. These 



were described theoretically at first by dedicated methods for the problem at hand, from 
Cohen-Tannoudji et al. (dressed atom picture) [13, 14] and Dalibard et al. (diagrammatic 
expansion) [15]. The extension of photo-detection in the spirit of Eberly and Wodkiewicz 
by considering two detectors with respective linewidths I\ and T 2 was impulsed by Knoll 
et al. [16] and Arnoldus and Nienhuis [17]. The expressions were of general validity, even 
though, due to their complexity, the authors still focused on the particular case of resonance 
fluorescence for illustration. The mathematical foundations, shaky in their initial develop- 
ment, were firmly established in the course of the following years [18-20]. The multiplicity 
of photons requires a careful time (T±) and normal (:) ordering of the operators [19, 20], 
and it was realized that it is the time ordering of (:Al J (Ti)A UJl (Ti)Al J (T2)A U2 (T2):) 

(2) 

which provides the physical two-photon spectrum S^ r (u)i,Ti;u)2,T2) = 

(7^^ (t'^a^ (t' 2 )}7 + [a(t' 3 )a(t' 4 )}) . Here, we have defined T + (resp. 7) to order the operators 
in a product with the latest time to the far left (resp. far right) [1]. Normalising this 
expression yields the sought time- and frequency- resolved two-photon correlation function 
9riT 2 ( w i> Tii w 2, T 2 ) = Sf^ 2 (ui, Ti; u> 2 , T 2 )/ [Sf^ (ui, Ti)SfJ (u 2 , T 2 )] ■ It is positive and finite, 
and reflects that frequency and time of emission cannot be both measured with arbitrary 
precision, in accordance with Heisenberg's uncertainty principle. The limiting behaviours of 

(2) 

<7pir 2 defined in this way are those expected on physical grounds: photons are uncorrelated 
at infinite delays, rim| T2 _r 1 |_ KX) <7p r {u)\ , T\ ', tu 2 , T 2 ) = 1 [21], and color-blind detectors recover 
the standard two-time correlators, limr 1) r 2 ^oofi , rir 2 ( Wl ' ^i! ^2,^2) = g^(Ti;T 2 ). Further 
generalisation to A^-photon correlations follows in this way, adding pairs of operators with 
their corresponding integrals [18, 22]. 

The actual computation of such j| ' r , however, have proved so far intractable for 
N > 2, even for simple single-mode systems, such as resonance fluorescence or the single 
mode laser [23]. The case iV = 2 is already demanding and thus some approximations were 
made to simplify the algebra [24, 25]. More recently, the resonance fluorescence problem was 
revisited without approximations but still for two photons and at zero time delay only [26]. 
The main reason for such limitations is that all the possible time orderings of the 27V-time 
correlator (X^a^i) . . . a^(t' N )]7+[a(t' N+1 ) . . . a(t 2N )}) result in (2N - l)!!2 Ar " 1 independent 
terms. Furthermore, each of these correlators requires the application of the quantum re- 
gression theorem 2N — 1 times. This growth of the complexity makes a direct computation 



hopeless for a quantity which is otherwise straightforward to measure experimentally, merely 
by detecting photon clicks as function of time and energy, a technology provided for instance 
by a streak camera [27]. 

In this letter, we present a theory of iV-photon correlations, that i) allows for arbitrary 
time delays and frequencies, ii) is applicable to any open quantum system and Hi) is both 
simple to implement and powerful. It consists in the introduction of N sensors to the 
dynamics of the open quantum system (noted Q in Fig. 1(a)). Each sensor of the set 
i — 1, . . . , N is a two- level system with annihilation operator q and transition frequency U{, 
that is matched to the frequency to be probed in the system. Its lifetime 1/Tj corresponds 
to the inverse detector linewidth. The coupling e^ to each sensor is small enough so that 
the dynamics of the system is unaltered by their presence, with (m) = (<? 4 q) <C 1. More 
precisely, calling 7q any transition rate within Q (either with internal or external degrees 
of freedom) linked to the field of interest a, the tunnelling rates e$ must be such that losses 



into the sensors and their back action are negligible, leading to e$ <C ^Y^q/2. Under 
this condition, we solve the full quantum dynamics of the system supplemented with the TV 
sensors. The latter play the role of the output fields A UJi (t) , but instead of formally solving the 
Heisenberg equations and expressing their correlations in terms of the system operators (as 
in the standard method exposed above), we compute directly intensity-intensity correlations 
between sensors, which is a considerably simpler task. The main result of this letter, which 
is demonstrated in the supplemental material, is: 

(at) ( „ „ v ,. (ni(Ti)...njv(Tjv)) 

9t 1 ..t n K u uTi\...\u3 N ,T n )= hm j^-r (1) 

1 N 6 U ...,e N ^o (m(7i)) . . . {n N (T N )) 

where the left hand side is the time- and frequency-resolved TV-photon correlation func- 
tion as defined previously [41]. The supplemental material establishes that, for open 
quantum systems described by Lindblad type master equations, (ni(Ti) . . . Un(Tn)) = 
p 1 " ' £ r N (2n) N Sl i r (o>i, 7\; . . . ; u^,T N ) to leading order in the e i: which proves Eq. (1). The 
equality is of general validity with no approximations or assumptions on the system. With 
this result, the complexity of computing g r r (oj\, 71; . . . ; un, Tn) is greatly reduced as no 
integral needs to be computed and the quantum regression theorem needs to be applied 
TV — 1 times only . For the important case of zero delay, g r r (ui, . . . ;u>n) reduces to 
a single-time averaged quantity. TV" degenerate sensors with frequency u and linewidth T 
also provide the iVth-order correlations of a single harmonic oscillator with frequency uj 



and linewidth T, corresponding to the case of correlations measured after the application 
of a single filter. This method is also useful to derive analytical results (as shown in the 
supplemental material). 

We now illustrate its efficiency and ease of use by applying it to the Jaynes-Cummings 
model [28], which is both an important and fundamental quantum description of light-matter 
interaction [29], is much more complex than resonance fluorescence as it also quantizes the 
light field [30] and is particularly suited to generate strongly correlated photons [31, 32]. 
Our method recovers exactly the known results for the Mollow triplet [24-26], and extends 
them effortlessly. 

At resonance between the light mode (a) and the two-level emitter (a) both with bare 
frequency u a , the Jaynes-Cummings Hamiltonian reads H = g(aV + aa^). The master 
equation that describes decay (7 a , 7 CT ) and incoherent pumping of the emitter (P CT ) has the 
form d t p = i[p, H] + [\L a + \L a + ^£ fft ](p), where L c (0) = (2cOJ - JcO - OJc) and 
p is the density matrix for the emitter/cavity system [33]. The new density matrix that 
includes the sensors, p sen , follows a modified master equation where the photonic tunnelling 
terms, H scn = X^il^W^ + e i{ a< H + a ^i)], are added to the original Hamiltonian, and the 
sensor decay terms J2 i=1 y^ ft (p S en) are added to the dissipative part. The level structure 
of the dressed states \n, ±) with n excitations is given by the dissipative Jaynes-Cummings 
ladder [33], which is shown in Fig. 1(b) at low pumping, P a = 7 CT , and in the strong- 
coupling regime with / ~f a < 7 a < 4g. This gives rise to the transition frequencies R^ = 
yng 2 — ( 7 °~ 7<T ) ± \J (n — \)g 2 — ( 7a ^ 7<7 ) between rungs for n > 2 with broadening 7„ = 
2(n — l)7 a + 7o- [33]. The Rabi splitting 2R, which arises from transitions |1±) — > |vac), is 
given by R = \J g 2 — ( 7a 4 7<7 ) with 7 X = (7 a + 7o-)/2. These transitions result in peaks in 
the power spectrum, as seen in Fig. 1(c) for the three cavity decay rates 7 a /5' = 0.01, 0.1 
and 0.5 that are chosen to correspond to cavities embedding superconducting qubits [34], 
atoms [35] and quantum dots [36], respectively. They all show the first rung transitions 
at ±R, the so-called Rabi doublet, and one can distinguish outer peaks at ±i?+ and inner 
peaks at ±-R^, up to the third rung for the best system (solid line) and to the second rung 
for the intermediate one (dashed line). In Fig. 1(d), we set the linewidth of the sensors T 
at a value around 72 and compute the two-photon correlation at zero delay, g^ (wi;^), 
between a photon with fixed frequency at the Rabi peak, 0J2 = R (solid arrow on the left of 
Fig. 1(b)), and a photon with variable frequency Ui which scans the spectral range (curved 



arrows). When the scanning frequency iO\ matches the second rung transitions that are 
precursors of the Rabi transition R, the probability of joint emission is enhanced relatively 
to other frequencies. The filtering then tracks photons in the cascades |2+) — > |1+) at 
Ui = R7, and |2— ) — > |1+) at —R 2 - This is a common feature to all three systems, which 
shows that even if broadening is too large to observe explicit features from higher rungs in 
the power spectrum, g r (ui, u) 2 ) allows to uncover them in the photon correlations. On the 
other hand, we obtain the expected strong suppression when the first photon is detected at 
the other branch of the Rabi doublet, u\ = —R. More features can be observed for the better 
systems such as dips at the two remaining transitions from the second rung, |2— ) — > |1— ) at 
Ui = —R 2 and |2+) — > |1— } at R%. In the best system, we can even resolve the dips for the 
third rung transitions at u\ = ±i? 3 • All these transitions do not form a consecutive cascade 
with the one we fixed and therefore have less probability to occur within the considered 

small time window I/72. 

(2) 
Instead of making a comprehensive analysis of g r specifics, we now turn to higher order 

correlation functions, such as the simultaneous three-photon correlations g r (oj\\ u 2 ; w 3 ), 

which are exceedingly hard to compute with previous methods. We fix two frequencies of 

detection at u 2 = R 2 and 003 = R (solid arrows on the right of Fig. 1(b)) and again let oo\ 

vary. A strong enhancement is also observed for all systems, now at ui\ = R% which monitors 

the cascade |3+) — > |2+) — > |1+) — > |vac) depicted in Fig. 1(b) and at oj\ = —R% which 

starts it with |3— } — > |2+). Other transitions show dips that are also clearly understood. 

This hints at the possible characterization of the level structure of an open quantum system. 

In general, however, one cannot draw conclusions from the zero-delay case only, in particular 

for small features, such as the small enhancement at u>i = —R 2 in g r (for the dashed line 

only) which is not necessarily a bunching peak and reveals itself in the r-dynamics to be 

antibunched, as discussed later. 

In Fig. 2(a), we explore another important aspect of g r , namely the dependence of 

correlations on the sensors linewidths, which is related to the complementary uncertainties 

in time and frequency. In the case r — > of perfect detectors, g Q — 1 for all iV with 

non-degenerate frequencies, since the complete indeterminacy in time leads to averaging 

photons from all possible time delays. For M degenerate frequencies out of N, photon 

indistinguishability results in M! ways for the sensors to measure the same configuration, 

that is, limp-s.ofl'r = -^- This limit has been misunderstood in the literature [42]. The 



effect has otherwise been reported for the case M = N = 2 by converting laser light 
into chaotic light with narrow filters [23]. The other limit r — > oo corresponds to the 
opposite situation of exact r-delay between photons of completely indeterminate frequencies. 
This is of more interest, in particular at zero time delay, which is the case of Fig. 2(a). 
For the Jaynes-Cummings system at low pumping, this recovers results derived by other 
approaches [37, 38]. 

The intermediate case of finite linewidth of the sensors is the most interesting. Features 
are the most marked when detector linewidths are of the order of those of the transitions 
involved, since the peaks of the spectrum are best filtered. Smaller linewidths (longer times) 
are to be favoured for bunching and larger linewidths (smaller times) for antibunching. One 
sees for instance in Fig. 2(a) that consecutive transitions, forming a cascade — such as those 
sketched in panel % (with three photons) or ii (with two photons) — show an enhancement. 
Conversely, the simultaneous emission from both Rabi peaks, in the configuration sketched 
as iv, is substantially suppressed, leading to strong antibunching. This observation with a 
microcavity containing a single quantum dot has been used to demonstrate the quantum 
nature of strong light-matter coupling [6] (with detuning to better separate the peaks). 
Further theoretical investigations with this formalism (to be discussed elsewhere) may allow 
to elucidate the nature of spectral triplets also observed in such experiments [6, 39, 40]. 

Figures 2(b-c) show an example of the r-dependence of the correlations, for the case 
r = 72, both at positive and negative delays. The configuration ii has the typical shape 
of a cascade between consecutive levels, with antibunching for r < 0, a step at r = and 
bunching for r > 0. This behaviour is well known, for instance from the biexciton-exciton 
cascade [9] . It is also observed for N photons in any consecutive transitions, such as is shown 
in % for three photons starting from the third rung. In contrast, the filtering of peaks which 
do not belong to the same cascade exhibit antibunching, as seen in iv for the two Rabi peaks 
or Hi for one of its three-photon counterparts: the order of the transition does not matter 
anyway and the cases ±r show qualitatively the same behaviour. These results are, to the 
best of our knowledge, the first computations of three-time frequency-resolved correlation 
functions. They are easily extended to higher orders (a fourth order example is given in the 
supplemental material). 

In conclusion, we have presented a theory to efficiently compute correlations between an 
arbitrary number of photons of any given frequencies and time delays. All three aspects of 



the detection, namely frequencies, time-delays and linewidths of the detectors, are needed 
to characterise meaningfully the system. The method allows to compute exactly, with low 
effort and for general open quantum systems, properties of output fields that are otherwise 
defined in terms of complicated integrals. Its ease of use enabled us to present the first 
computation of three and four time-resolved and frequency-filtered correlation functions. 
Its application will allow the interpretation of experiments which are routinely implemented 
in the laboratory but which lacked hitherto an adequate and tractable theoretical support, 
and to design new ways to unravel and/or engineer the quantum dynamics of open systems. 
EdV acknowledges support from the Alexander von Humboldt foundation; AGT from the 
FPU program AP2008-00101 (MICINN); FPL from the Marie Curie IEF 'SQOD' and the 
RyC program; CT from MAT2011-22997 (MINECO) and S-2009/ESP-1503 (CAM); MJH 
from the Emmy Noether project HA 5593/1-1 and from CRC 631 (DFG). 
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FIG. 1: (Color online) (a) Scheme of our proposal to compute iV-photon correlations between 
photons emitted at different times and frequencies from an open quantum open system Q. N two- 
level systems of ascribed frequencies are weakly coupled to Q and serve as correlation sensors at 
these frequencies, with their decay rate providing the detector linewidth. (b) Dissipative Jaynes- 
Cummings ladder up to the third rung with two of the cascades probed in panels (d) [with two 
sensors] and (e) [with three sensors]. Solid arrows show the fixed frequencies. Curved arrows 
show the scanning frequency uj\, at the transitions where the joint emission is strongly enhanced 
(dashed) or, on the other hand, suppressed (dotted), (c) Power spectra of emission probed by weak 
incoherent excitation (P ff = J& = O.Olg) for three cavities of decreasing quality 7 a = 0.01 (solid), 
0.1 (dashed) and 0.5<? (dotted), (d) Two- and (e) three-photon correlations at zero delay for the 
three cavities, with sensor linewidths T = 72 (solid) and 72/2 (dashed and dotted). 
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FIG. 2: (Color online) (a) Two- and three- photon correlations at zero delay as a function of the 
sensor linewidth V, with frequencies of detection as shown in the insets i-iv. (b-c) r-dynamics of 
the correlation functions with T = 72 for, (b), two photons in the configurations of insets ii and iv 
and, (c), three photons in the configurations i and in. Positive r corresponds to detection in the 
order from top to bottom of the ladder. Parameters: P a = 7^ = O.OI5, 7 a = O.lg. 
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We 



prove 

9r 1 '..r N ( UJ ^ T ' 



the equality between the experimentally motivated correlation function 
i;...;un,Tn) — defined as 2iV-time integrals — and the intensity correlation 
(ni(Ti) . . . riiv(Tjv)}/((ni(Ti)} . . . (niv(Tiv)}) between sensors coupled to the system, to lead- 
ing order in the coupling. We also further illustrate the method by calculating correlations in the 
Jaynes-Cummings model up to fourth order. 



PROOF OF THE EQUIVALENCE BETWEEN 
THE SENSING AND THE INTEGRAL METHODS 



Let us assume a quantum system described by a set 
of operators a, a, etc., acting in a Hilbert space H. In 
second quantization, these operators define annihilation 
operators in the Heisenberg picture. The system can 
be fully Bosonic, Fermionic or a mixture involving any 
number of operators. All single-time quantities can be 
obtained from correlators of the type (a) ^ a? tj^' n a 9 . . .) 
with [a, v, rj, 6, etc., integers. Let us call O the set of 
operators the averages of which correspond to the corre- 
lators required to describe the system, i.e., O includes all 
the sought observables as well as operators which couple 
to them through the equations of motion. In the fol- 
lowing, we assume, without loss of generality, that a is 
the mode of interest, the correlations of which are to be 
computed in time and frequency. 

We prove the case N = 1 first, which corresponds to 
the power spectrum, then N = 2, which corresponds to 
the most important correlation function. The proof ad- 
mits a straightforward generalization to higher N. The 
proof proceeds by computing separately the integral ex- 
pressions on the one hand and the intensity correlations 
between sensors on the other hand, and showing that 
they are equal to leading order in the couplings to the 
sensors. We assume the steady state case for simplicity, 
with little loss of generality. 



N — 1, power spectrum 



Integral method 



The single-photon physical spectrum intro- 
duced in the main text as <Sp (wi,Ti) = -^ x 

can, through convolutions, be put in the form of an 



uncertainty in the time of detection [1] : 



SgVi.Ti) 

where 

4 1 1 ) (wi,*i) = -» 



Ti 



dtie- r ^ T -^)4V(«i.*i) (!) 



dT ie -^ Tl e- lu,lTl 



(^(iiMii-n)) 

(2) 
contains the uncertainty in the frequency of detection [2]: 

The 



£Wi) = /! 



DC' 

oo 



.(1 



dwiE^K.tiJi- 



(iJ^ + K-u,!) 2 

kernel of this expression corresponds to the case of a 
perfect detector, Ti = 0, known as the Page-Lampard 
quasi-spectrum of emission Sq (wi, ii) [3]. The results of 
Ebcrly and Wodkiewicz [1] show that the time-dependent 
physical spectrum (1) is i) always positive, whereas 

Sq (wi,ti) is not in general, and it) finite, even in the 

_ V M 



= Sp (wi,ti), whereas 



steady state Sp ' (wi , T\ 
S^\uji,Ti — >• oo) diverges. 

To compute Eq. (1), we only need to obtain the two- 
time correlator (a' (t\)a(ti — Ti)). For any two operators 
X and Y acting on H, we define the vector Vx,y(t) as: 



v x,y(t) 



/ (A(o)y(o)) \ 

(X(0)a(r)F(0)> 

(Jf(0)ot(r)y(0)> 

(X(0)(«ta)(r)y(0)) 



V 



(3) 



/ 



where X and Y, in the steady state, sandwich the oper- 
ators of O taken in some order, which will be kept for 
the remainder of the text as starting with the sequence 
O = {1, a, a' , a^a, . . .}. 

From the quantum regression theorem, one can define 
for O a matrix M which rules the dynamical evolution 
of v X y: 

d T v x ,Y(r) = Mv x ,y(t) , (4) 

with solution Vx,y(t) = e Mr vx.y(0). The steady state 



of the system is then fully given by: 



v ss = lim vi.i(r) = lim e 



Mt I 



(5) 



since O contains all the relevant observables of the sys- 
tem. Here we have chosen the vacuum as the initial con- 
dition. Since we employ the standard assumption of a 
unique steady state, the initial state does not matter and 
all the information is encoded in e Mr . 

We now define two matrices, T±, which, when acting 
on vx,y(t), introduce an extra a) for T + and an a for 
T_ between X and Y, keeping normal ordering: 



T+Vx,y(t) = 



and 



T_v x ,y(t) = 



(X(0)«t(r)F(0)) \ 
(X(0)(«ta)(r)y(0)) 

(X(0)a^(r)Y(0)) 
(X(0)(a^a)(r)Y(0)) 

V '■■ J 



(X(0)a(r)Y(0)) 

(X(0)a 2 (r)y(0)) 
(X(0)(ata)(r)y(0)) 

(X(0)(ata 2 )(r)r(0)> 



(6) 



V 



(7) 



/ 



These matrices always exist, in infinite or in truncated 
Hilbcrt spaces (where, if truncation is to order n, a n is an 
operator in O but a n+1 = 0). For instance, if the mode a 
is a two-level system, the vector vx.y(t) consists of the 
first four entries in Eq. (3) only, since (V^a v = if /i or 
v > 1. Then, these matrices read 



T + = 



1 


V o o o o / 



and T_ = 




1 

V o o o o / 



(8) 



With these definitions, the correlator (aJ(ti)a(ti — t\)) 
with t\ > is the first element of T + Vi j0 (ri): 



(at(t 1 )a(ti-r 1 )) = [T +e M ^T_v ss ] 1 , 



(9) 



where we have used [•••]* to denote the zth element of 
a vector. The power spectrum in its integral form is 
therefore given by: 



,(i) 



1 



S^(wi) = -K 



T + 



-1 



'M + (-iwi-^)l 



-T_v s 



(10) 



with 1 the identity matrix. 

Sensing method 

We now consider two sensors q, i = 1,2 with 
linewidths I\ coupled to the system with strength Si such 
that the dynamics of the system is probed but is other- 
wise left unperturbed. This requires the tunnelling rates 
Si to fulfil two conditions: the losses into the sensors must 
be negligible, ief/Ti <c -jq and so must be the back ac- 
tion of the sensors into the system, 4e 2 /7g <C 1^, where 
7q is the smallest system decay rate. These conditions 
both lead to e, <C ^JTc/q] 7 !. We then introduce a sens- 
ing vector w of steady state correlators, by multiplying 
q tMi^i ( -tP2^ 2 with the opcratoi . s in q- 



( 



w[/ii^i,^ 2 ^2] = 



V>1 M ^2 <>2 



) \ 



fci <*'&'<%' a) 



\ 



(11) 



/ 



where the indices /Zj and Vi take the values or 1. In the 
regime under consideration, the population (?/ft) <C 1 
and the equations of motion are valid to leading order 
in £ 12 : 



Ti r 2 

(9 t w[/xi^i,^ 2 ^2] = {Af + [(/ii - !yi)iwi - (^i +^i)— + (^2 - ^)iw 2 - (/-*2 + t'2)-^-]l}w[ J uii/i,^ 2 z/ 2 ] 

+ ^i(i£iT+)w[0z/i, / u 2 f 2 ] +z/i(-ieiT_)w[yUiO,,u2^2] + /U2(«£2T , +)w[ / ui^i,0f 2 ] + f 2 (-ie 2 7 1 _)w[)Uii/i,^ 2 0] , (12) 
and can be solved recursively: 

-1 



w[/ii^i, ^2^2] = i 

M + [(Mi - ^i)iwi - (Mi + v i)^t + (M2 - ^2)^2 - (M2 + t / 2)-#]l 

x | y ui(zeiT + )w[0^i,^2^2] +^i(-i£iT_)w[yUiO, ^2^2] + M2(^2T , +)w[ / ui^i,0za 2 ] + v 2 { -is 2 T-)w [^1^1, M2O] j • (13) 

Higher order terms will cancel exactly in the vanishing coupling we will assume later and thus do not need to be 
included here. Besides, unlike the leading order term, the higher order ones depend on the modelling of the sensors 
(as two-level systems, harmonic oscillators, etc.) and on the system itself. 



The spectrum of emission of a is given by the average 
population, in the steady state, of any one of the two 
sensors, say, the first one: (m) = (TiCi)- Its equation 
of motion reads d t (ni) = —Ti(ni) + 2-ft(i£i(9 L a' f )), and 
with the above notations, is therefore given in the steady 
state by: 



J- 1 



ieiT+w[01, 



(14) 



Using the solution Eq. (13), the correlator of interest for 
the spectrum reads: 



w[01,0,0] = 



-1 



M+[-twi-^]l 



H £l T_)v s 



(15) 



identity: 



(m> 



2 A 



1? 



r. 



M ■ 



-t«i - ^]1 



-T_v s 



Ti 



I (2^)^ 1 , (wi). (16) 



N = 2, two-photon correlations 



Integral method 



Equality of the integral and sensing methods 

The proof is now complete since, to leading order, we 
find that Eq. (10) and Eqs. (14-15) provide the claimed 



The case N — 2 brings with the multiplicity of photons 
the conceptual difficulty of time- and normal-ordering. 
It was discussed in the text that the proper definition 
yielding a physical two-photon spectrum reads [4, 5]: 



rir 2 



Ti 



x e^^-^h^^-^^T-ia^t^a^t'^T+Ht'^a^)}) . (17) 



In analogy with the case N = 1, it can be put in the form 



S^r^wt^u^Ti) = 1^2 J' ' dti f ' dt 2 e- r ^ Tl - t ^e- r ^ T2 - t ^ ) iT2 {u 1 ,t 1 ;u2,h), (18) 

J — oo J —oo 

isolating the two-photon quasi- distribution: 



.(2) 



21? 



ri r 2 ( w i^i; w 2,*2) - ^y 



dr\dr 2 e ? Tl e 2 



x e-^le^HT-laHti - ^ {t 2 )]T+[a(t 2 - T 2 )a{h)]) + e"*" ^ (T-[<J(h)a\t 2 )]T + [a(t 2 - r 2 )a(h - n)]>] , (19) 



which, like the quasi-spectrum, can be negative and is thus not a physical spectrum. 



To proceed with the calculation, let us separate the r = 
T 2 —Ti two-photon correlation function between its r = 
and r > terms: 



Sr-Vr^wi;^,"!-) 
r 2 r c (2) 



(2) 



e- l * T Sg ra (wuW2) + AS^ r2 (wi;w2,r) , (20) 



with 



Sp ir2 (wi;w 2 ) = Tir 2 / cft 2 / dti 

J —OO J— OO 

x e - r 1 ^-* 1 ) e -r 2 (T 1 -t 2 ) s (2) r2(u;i)ii;w2)i2) + [1 ^ 2] 



(21) 



and 



r r 2 



r'A 



A4 2 1 ) r 2 (wi;w 2 ,r) = r 1 r 2 [ ' dh [ * dti 

JTx J -co 

x e -rim-t 1 ) e -r a (T 2 -t 2)s g) ra(wijti . W2it2)i (22) 

where [1 -H- 2] means the interchange of sensors 1 and 2, 
that is, permuting uj\ -H- u; 2 and Ti <->• r 2 . 

To compute these quantities, it is enough to consider 
Sp ' r (wi, ti; w 2 ,i 2 ) for i = i 2 — ii > since the inverse 
order is given by the exchange [1 -H- 2]. Therefore, we 
restrict the integration to ordering of the time variables 
where t\ — t\ <t\ < t 2 . The fourth variable yields three 
different domains of integration: 



(l) t 2 -r 2 <ti-n<h <t 2 , 

(2) tl -Tl < *2-T 2 <il < t 2 , 

(3) tl -Tl < tl < t 2 -T2 < t 2 . 

For each of them, there are two different correlators ap- 
pearing in £( 2 ); one with the factor e _lW2T2 e lWlTl , the 
other with e - *'*' 2T2 e -w *' lTl . They will be respectively re- 
ferred to as C/ ia \ and Cu^, with i = 1,2,3 depending 
on their domains of integration. This gives rise to six 
integrals which we shall denote T-i ia \ and lub) ■ 

From this discussion, we can find a general expression 
for the complexity of the integration method in terms 
of the various domains of integration and the different 
correlators to be considered. The number of independent 
time ordering is (27V— 1)!! and the number of independent 
terms in EW is 2^/2 (we divide by 2 because half are 
complex conjugates of the other half). The total number 
of independent time integrals and correlators is therefore 
{2N-l)l\2 N - 1 . 

The first correlator we need, Cn a ) = ( a ^(ti — 
Ti)er (t2)a(ti)a(t2 — T2)), is the first element of the vec- 
tor r+Vxi,Yi(t) with X\ = a) [t\ — t\) and Y\ = 
a (t!)a(t 2 -V 2 ). We obtain w Xl , Yl {t) = e Mt v XliYl (0) - 



e M 'T_v XliF2 (Ti) with Y 2 



a(h 



T2 



In turn, 



Vx u Yi(n) = e M ^v Xl ,y 2 (0) = e MTl T + Vi,r 2 (t') is ob- 
tained with Y 2 = a(t 2 — t 2 ) and t' = t 2 — t\ — t. Finally, 



we get Vi iY2 (t') 
everything together, we get: 



e Mt 'vi,vM 



,Mt' 



T_v s 



C(la) 

(16) 



r,e M *T T e MTl T+e M *'T_v ss 



Putting 



(23) 



with correspondence between upper and lower indices 
with the sign. Repeating this procedure for the other 
domains of integration, we also get: 



C(2a) 
(26) 



rr, Mtrp -Mt Mr 2 rr Mt" rp s 

l + e lzre e 2 i„e l±v 



(24) 



where we defined t" — t + t\ — t 2 (going from to 00), 
and 



C(s.) = [T + e MT2 r_e M *e- MT2 T T e MTl T ± v ss ]. 



(25) 



Integral method at r — 



(a) and (6) read: 

rir 2 1 



I(la) — 



(i6) Ti + T 2 (2tt) 



T, 



-1 



M + {-iuj 2 -T 1 - v -f)l 



xT T 



-1 



^)1 



M + (±icu 1 - iuj 2 r,+J 

x T± — = T_v s 

M +(-»W2-^)l 



(26) 



The second correlators Ci 2i \ lead to: 



r^ 1 

(26) Ti + T 2 (2tt) 2 



T, 



-1 



xT x 



'M + (-iw2-ri-^)l 

-1 



M + (±iui - iuj 2 - ti f ia )l 
x T_ — = T+v s 



M + (±iwi - ^)1 
And the third correlators C/3^ lead to: 



(27) 



rir 2 1 



•<36)' ri + r 2 (2^) 2 

x T_ 
xT 5 



T, 



-1 



M + (-iw 2 - Ti - ±f)l 
-1 



M-Til 
-1 



M+(±iwi-^)l 



-T+v s 



(28) 



The total correlation function follows from twice the 
real part of the six previous integrals summed over and 
exchanging photons: 



j=l,2,3 



) i-%6) 



+ [1 ++ 2] . (29) 



Integral method at r > 



:(2) 



The finite time-delay contribution ASp r (u>i ;uj 2 ,t) 
requires different domains of integration only for the vari- 
ables t 2 , now ranging from 1\ to T 2 , and t = t 2 — t\ now 
ranging from t 2 — T± to 00. As a result, the integrals in 
Eq. (22) depend on r. The integrals on the correlators 
C {ia) and C{2a) , that we note AX(i a ) and AZ( 2 a> give sim- 

(16) (26) _ (16) (26) 

ilar results as the corresponding I{ ia ) , but they acquire 

(i6) 

the T-dependence in the form of a factor (ri + T 2 )F(t), 
with 



We now turn to the zero time delay contribution 



j(2) 



(uji;uj 2 ), which, according to Eq. (21), is given by 
integrating the correlators (Eqs. (23-25)) over their cor- 
responding domains, changing variables as needed. For 
instance, the integrals of correlators Cm) require the 
change of variables t\ —¥ t and t 2 — > t' (both extending 
from to 00). The final expressions for the two integrals 



[M+(-iw 2 + i 2 2 -)l]r _ I 

J-(r) = e- r2T = , (30) 

that is to be inserted in Eqs. (26), (27) after the first 
matrix T + . The integrals on C(3a>, on the other hand, 

(36) 

arc not so straightforward. They are to be separated 
into two parts: one where t 2 — t 2 < T\ 1 the other 



one t 2 ~ r 2 > T\. The first part, with integrals 
J T 2 dt 2 J T dt J t _ T dr2 J dT\{. . .), gives rise to a quan- 
tity similar to AZ<;a) (r) with i = 1,2, in that its r- 

<ib) 

dependence also consists in the factor (Ti + T 2 )F(t) in- 
serted after the first matrix T + in Eq. (28). For this 
reason we note it AI< 3 a) (t). The second part, with inte- 



grals j T 2 dt 2 j t _ T dt j Q 2 1 dr 2 L dr\{. ..), yields two 
more contributions: 



AX (3a) (r) = ^| 

(3/3) V^Tt) 



T+Z(r) 



-T±v s 



(31) 



M-Til + M + (±iwi - 
where we introduced the r-dependent matrix 

Z(t) = dtj dr 2 e- r ^^ 

Jt x Jo 

X e ( M +(-^2- ! #)l)r2 r _ e M(t 2 -T 1 -T 2 ) _ / 32 \ 

Z(t) can be calculated for each element T^} <G {0, 1} of 
the matrix T_ : 



^A T ) 



771 TTi—lrrikl TP T?— 1 

2 



2T E 

p,k,l,q P 1 * 



^(m p -iu 2 + ^-)T _ y e {m q +T 2 )T _ -y 



m p - ico 2 + 



}, (33) 



where E is the matrix of eigenvectors of M, that diago- 
nalises it: M^ = ^ E ip m p E~^, with m p the eigenval- 
ues. 

Gathering terms with the same r dependence defines 

^ I ( T ) = Z)j=i,2,3 A2: (ia)(T") + AI( j6 )(t) which enters 
in the final result: 

A4?r a («i;«2,r) = 23?[az(t)+AZ (3q) (t)+AZ (3W (t)" . 

(34) 



This solution relics on w[ll,01] which can be expressed 
in terms of three lower order correlators: 



-1 



w[ll,01] = 

M + (-»w2-ri-^)i 



x |-ie 2 T_w[ll, 00]-ieiT_w[10, 01]-MeiT+w[01, 01]} 



(37) 



each of which is given by: 
w[ll,00] = 



M-ra 

x |i £l T+w[01, 00] - kiT_w[10, 00] } , (38) 



410,01] 



•1 



M+(iui -iuj 2 - El ^)l 



x I - i£ 2 T_ w[10, 00] + i£iT+w[00, 01] } , (39) 
and 

w[01,01] = 



-1 



M + (—iu>i — iuj 2 



ri+r 2 



x I - i ei T_w[00, 01] - i£ 2 T_w[01, 00] } . (40) 

Finally, we apply a last time the recursive relation 
Eq. (13) to find the three different correlators involved 
in the previous expressions: 



w[10,00] = 
w[00,01] = 
w[01,00l = 



-1 



-ie x T + w ss , 



(41a) 



M + (iwi-^)l 

, f + , 7 1 r H £2 T_)v B8 , (4ib) 

F - T -(- J£l T_)v ss . (41c) 



M + (-iwi-^)l 



Sensing method at r > 



Sensing method at r = 

The intensity correlations between two sensors, 
{n\n 2 ) — (^^1^2^2)1 have the equation of motion: 



d t (nin 2 ) = -(Ti + r 2 )(nin 2 ) 

+ 23? 



^(qci&a*) +i£i(<ri4^2a t ) 



(35) 



This leads to the steady state solution 
2 



(nin 2 ) 



ri + r 2 



-3? 



i£ 2 T+w[ll,01] + [!«->■ 2]. (36) 



We now consider the case where the second photon is 
absorbed by sensor 2 with some delay r > after a first 
photon is absorbed by sensor 1. The correlator of interest 
is (Hi(0)n 2 (r)), with equation of motion: 

9 T (ni(0)n 2 (r)) = -r 2 (ni(0)n 2 (r)) 



23? 



^MOXWXr)) , (42) 



with the initial condition in the steady state 
(ni(0)n 2 (0)) = (nin 2 ). This solution relies on 
(ni(0)(c 2 aT)(r)). To compute it, we introduce a vector 
analogous to Eq. (43) but now consisting of two-time 



correlators: 



following order: 



w'[ll,/i 2 ^2]( r ) = 



( ni (0)(4^a)(T)) 
(m(0)(^^ o t)(r)) 
MOX^^otoXr)) 



V 



(43) 



/ 



With this definition, (ni(0)(? 2 er)(r)) is the first ele- 
ment of the vector T + w'[ll,01](r). The r-equation for 
w'[ll, 01] (r) reads: 

d r w'[ll, 01](r) = [M + (-1W2 - y )l] w'[ll, 01](r) 

- i£ 2 T_ w' [1 1, 00] (r), (44) 

with w'[ll,00](r) = e Mr w[ll,00] with initial condition 
w'[ll,01](0) = w[ll,01] in the steady state. After some 
algebra, one arrives to the solution: 

w'[ll,01](r) = e[ M+( - 4W2 -^ ):l ] r w[ll,0i] 

-(-te 2 )y(T)w[ll,00], (45) 

in terms of a matrix y(r) defined clcmcntwise as: 



ya(r)= E 



E R~ l T kl E, Fr 1 



p.k,l,q y y 



U0 2 



e (ra p -!u 2 -^)r _ e r; 



* T } . (46) 



Substituting this expression into Eq. (42) and solving it, 
we obtain: 



(ni(0)n 2 (r))=e- l2r (n 1 n 2 ) 



23? 



i£ 2 T+J"(r)w[ll,01] 



23? 



elr + Z(r)w[ll,00] , (47) 



where the matrices F{t) and Z(r) are those introduced 
in the previous section, namely, Eqs. (30) and (33), re- 
spectively. 



Equality of the integral and sensing methods 

We complete the proof by showing that the results from 
the integration and the sensing methods are the same to 
leading order in the coupling s. 

First, the case r = 0. The final expression for (nin 2 ) 
is obtained by inserting the solutions for the correla- 
tors (38-40) into Eq. (36). This leads to the same results 
as Eq. (29), with the integrals appearing precisely in the 



_2„2 



{n 1 n 2 ) 



TiT 



23?{: 



11 2 



-(27T) 



•(3b) 



'-(3a) 



1, 



[1++2] 



(2a) 

2 2 
£ 1 £ 2 

rir 2 



-i, 



(la) 
x2c(2) 



-x, 



(16) 



'-(2b) 



(27r)^ r > i;W2 ). (48) 



Second, the case r > 0. For ease of comparison, we 
rewrite the term AZ in term of the vector w[ll, 01] as: 



Al(r) = 



(2ny 



T + T(i 



1 



e\{-ie 2 ) 



,[11,01] 



(49) 



It is then clear that this expression is equal, up to a 
constant factor, to the second line in the expression for 
(ni(0)n 2 (r)), Eq. (47). Similarly, the term AI( 3a )(r) + 
A1( 3 ^)(t) in Eq. (31) can be rewritten as: 



£Z {3a) ( T ) + Al m ( T ) 



rir 2 



(27T) 



T+Z(r)lw[ll 



(50) 

and related to the third line in Eq. (47). All together, 
we can therefore conclude that, to leading order in the 
couplings: 



<ni(0)n 2 (r)> = ^$-(2n) 2 S^ ra ( Ul]U2 ,r). 

t it 2 



Final remarks 



(51) 



This proof can be generalised to A-photon correlations 
and/or for finite Ti-time dynamics (instead of a steady 
state) by repeating these procedures linearly in the num- 
ber of sensors and integrals. There is no conceptual dif- 
ference brought by the higher number of variables, but 
notations become heavy and for the sake of clarity, we 
have illustrated the proof in the simplest, as well as most 
relevant cases, of N = 1 and 2. Also, nothing in the proof 
relies on the choice of sensors as two-level systems, which 
has been made for convenience. As we always examine 
crossed correlations between them, they could also be, 
e.g., harmonic oscillators, and provide identical results. 

Together with Eqs. (36-41), Eq. (47) provides a semi- 
analytical result that can be used directly for computa- 
tions. Although the Hilbert space is not enlarged when 
using these formulas, they are however awkward to use 
and set up. Also, the growth in the number of correla- 
tors has the same power dependence on the maximum 
number of excitations allowed in the system than when 
including the sensors explicitly (it is linear in the Jaynes- 
Cummings model) . The number of correlators increases 
like 4^ when including N sensors. Benchmarks for N = 2 
show that the many matrix operations (inversions and 
multiplications) involved to evaluate the formulas are 



more costly than solving linear equations as required 
when including explicitly the sensors. Although this is for 
a larger set of correlators in the latter case, optimisations 
such as LU decomposition make sensors a more efficient 
as well as a conceptually simpler approach. If using the 
semi-analytical formulas turns out to be more effective in 
a particular context or for larger TV, similar results can 
be derived for (ni(0)ri2(ri) . . . tijv(tjv-i)) by generalizing 
Eq. (13) with TV sensors to obtain w[/xi^i, . . . , hn^n] re- 
cursively. 

Finally, the case of TV identical sensors reproduces 
exactly the TV-photon correlations, g^ N \ from a sin- 
gle harmonic sensor (full correlations of the output of 
a single filter). This can be shown by comparing the 
presented derivation with TV two-level sensors with one 
where the system is coupled to a single bosonic sensor 
with associated w vectors of the type w[n,m] (where n, 
m = 0,...N). The results are also seen to be identi- 
cal to those obtained by substituting u>i = 102 = w and 
r"i = T2 = r in the formula for <?p r (u>i; u>2,t). 



FURTHER APPLICATION TO THE 

JAYNES-CUMMINGS MODEL 

The sensing method was illustrated in the text up to 
three frequencies and for various time delays. Here we 
provide a supplemental example up to four-photon cor- 
relations. 

Such high-order correlations are not intuitive to visu- 
alise in their most general representation, given that they 
convey more information and of a much deeper character 
than single-photon observables. Photons are emitted at 
all energies and some correlations for particular energies 
other than ±i?„ are suppressed or, on the contrary, en- 
hanced, meaning that more complicated processes than 
simple relaxation take place. We reserve to future works 
the presentation of how such new processes of emission 
can be identified in the study of frequency resolved cor- 
relations, already at the two-photon level, and how these 
may find new applications to optimise quantum emitters. 
Here, to keep the discussion succinct, we will focus on 
the most important processes only, where TV photons are 
detected at precisely the Jaynes-Cummings transitions, 
that is, we disregard the correlations where one or more 
photons have an energy which does not correspond to a 
transition in the ladder. 

In Fig. I, we compare two, three and four-photon cor- 
relations of photons with energies corresponding to the 
possible placements of detectors over all possible tran- 
sitions. There are 2 2N ~ 1 configurations, half of them 
being symmetric with the other half by the interchange 
of upper and lower polaritons in all rungs. We need only 
consider, therefore, 2 2N ~ 2 cases, which are displayed in 
Fig. I representing only the half which detects the upper 
polariton in the first rung. For instance the leftmost case 



in panel (c) corresponds to setting four detectors at the 
energies wi, . . . , W4 probing the transitions |4+) — > |3— ), 
1 3+) -¥ |2-), 1 2+) -> |1-) and |1+) -)• |vac). As in this 
sequence of detection, polaritons have to swap branch in 
all rungs, the emission is unlikely and the corresponding 
coincidence is strongly suppressed. 



As discussed in the main text, finite r are important 
since correlations may be maximised at nonzero time- 
delays, when the dynamics of relaxation synchronises 
with detection. It is however difficult to find the opti- 
mising values for TV — 1 independent degrees of freedom 
when measuring TVth order correlations. We show here 
that the simplest approximation to fix all delays at zero 
already leads to useful results which contain the gist of 
the dynamics. An absolute value of photon correlations 
has little meaning in itself. It is when compared to other 
correlations in alternative configurations that a physi- 
cal meaning can be identified and quantified. Figure I 
shows how, even at equal times, the detection of photons 
at energies that correspond to a cascade of the Jaynes- 
Cummings ladder results in giant bunching. These are 
the points on the right of each panel. On the opposite, as 
previously described, when the detectors are arranged to 
click in the sequence that least correspond to a cascade, 
that is, alternating the type of polariton each time the 
system goes one rung down, a corresponding giant sup- 
pression is obtained. This is an actual antibunching in 
the case of two-photon detection (a) , while with a higher 
number of photons, the values obtained are larger than 
one but, again, when compared to the relative values of 
other transitions, reveal a giant suppression of correla- 
tions of over five and ten orders of magnitudes in three 
and four-photon counting, respectively. 



Another remarkable behaviour of these figures is the 
emergence of a classical behaviour with the increasing 
number of detected photons, powered by combinatorial 
growth. While correlations are markedly distinct and 
varying abruptly in the extreme, low-entropy situations 
(on both sides of the horizontal axes) , the large number 
of intermediate configurations smoothes out the quan- 
tized character and yields a gradual and milder variation 
as one quantum in the chain of detections is shifted from 
its precise expected value. The larger the number of pho- 
tons, the faster is this transition from a discrete, stair- 
case behaviour to a smooth continuous one. These are 
the transitions shown in black (also with smaller arrows 
in (c)). These results also reveal that proper sequencing 
of the detection allows to isolate and magnify its quan- 
tum character, even when dealing with a large number of 
photons. This only hints at the rich physics unravelled 
by TV-photon correlations and at the applications they 
could bring about. 



(a) 



|2+> 
12-5 

|i-5 



|vac) 











3 




IT V 


\ i 






CM 

\ ° u ir \r ir 




(b) 



2 3 

Configurations 



|3+) 
|3-) 

|2+> 
|2-> 

|l+> 

| vac) 







t 








1— 1 

3 




I 


^ r ^ f m ' 


f n • ▼ ^ 


i r 


11 " . T 


T i f ir 


w 


. ▼ " ir 


CM 

3 












CO 

3 , 















CO 

3 

(N 

3 



3 






iff 




' 


.» • 


1U 

io 5 

10 4 






» « 


3 

10 

2 

10 
10 


..- •■• 


• *•■ 


■■• 




1 4 7 


10 


13 16 



Configurations 



(c) 



3 

co 

3 

CM 

3 



|4- 



2+ 
2- 

1 + 

1- 



|vac) 




10 b 
10' 
10 f 



1 10 

?r io 



-tt*- 



.•■•■•*' 



1 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46 49 52 55 58 61 64 

Configurations 

FIG. 1: (Color online) (a) Two-, (b) three- and (c) four-photon correlations in the Jaynes-Cummings model with detection 
at the frequencies corresponding to transitions in the ladder. For visibility, whenever the arrows would overlap, a small shift 
has been introduced to assist the eye in tracking the starting and ending points. By probing transitions in a cascade or, on 
the contrary, from different de-excitation routes, the correlations between the detected photons vary over 3, 5 and 10 orders 
of magnitudes at the two, three and four-sensor level, respectively. Note that the vertical axes are in log-scale, so even if the 
variations appear moderate, they are locally important. Only transitions detecting the upper polariton |1+) — > |vac) in the first 
rung are shown, those detecting the lower polariton |1— } — > |vac) are reconstructed by swapping upper and lower polaritons in 
all rungs. Parameters: 7 a = 7 CT = O.OOlg, r = 72 = 0.003(? in the limit of vanishing incoherent pumping of the emitter, P CT — > 0. 
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